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Abstract: 

We use the BNL E821 measurement of — 2, the anomalous magnetic moment of the muon, 
to normalize, within a supersymmetric GUT framework, constrained MSSM (CMSSM) pre- 
dictions for processes that violate charged-lepton flavour conservation, including fi — > 67, 
— > e conversion and fi^e^. We illustrate our analysis with two examples of lepton 

mass matrix textures motivated by data on neutrino oscillations. We find that /i — > 67 may 
well occur at a rate within one or two (two or three) orders of magnitude of the present 
experimental upper limit ii — 2 is within the one- (two-)standard deviation range indicated 
by E821. We also find that /x — > e conversion is likely to occur at rate measurable by MECO, 
and there is a chance that — > /i^e^ may be observable in an experiment using an intense 
proton source. 
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1 Theoretical Framework 



The observation of neutrino oscillations [jl], Q implies that the individual lepton numbers 
Le,fi,T are violated, suggesting the appearance of processes that violate charged-lepton- 
number, such as /i — 67, /x ^ 3e, /i — ^ e conversion on heavy nuclei 0, r — /i7 
and Kl fie [H. The present experimental upper limits on these processes are 
B{fi 67) < 1.2 X 10-^1 §, ^ e+e+e") < 1.0 x IQ-^^ ^(^-Tz ^ e'Ti) < 

6.1x10-13 |]^^(^^^^^) < 1.1x10-6 g and ^ /i±e^) < 4.7x lO-^^ |g_ On the 

other hand, in minimal GUT models where the small neutrino masses are generated 
by the see-saw mechanism with massive singlet neutrinos un, and there are no new 
lighter particles, the amplitudes for charged-lepton-flavour violation are proportional 
to inverse powers of the heavy singlet neutrino mass M^^, and the rates for rare decays 
are extremely suppressed |1TT| . 



However, the observation of an apparent discrepancy between the measured value of 
the anomalous magnetic moment of the muon, = {g^ — 2)/2, and the value predicted 



in the Standard Model |T2[ suggests the appearance of new physics at the TeV scale 
in the lepton sector, with superseymmetry being one of the favoured options |T^, 0]. 
Moreover, there is a striking resemblance between the effective operators that generate 
fi e'j and 5a^: 

J^eff = e^£,a^,F'^^(A^f + Aff P^)£, (1) 

resulting in 



W 



Br(/.-^e7) = ^(|Ai^^^P + K/^n (2) 



and 



Sa, = ^iAir + Ar). (3) 

Hence, if these quantities are dominated by either the A^j or the A^, there is a direct 
relation between them: 



B{fx 67) = X {5a^f X e^ (4) 



where the lepton mixing factor e = A^^'^^/A^/^'^^. We see explicitly from (^) that 
the apparent measurement of 6a^ enables the rate for — >• 67 to be predicted, in 
the context of any model of lepton flavour violation motivated by the observations of 
neutrino oscillations which is able to predict e. 



One example of a theory where this connection can be made is supersymmetry KB 



1^. In a supersymmetric model, the amplitudes for processes violating charged lepton 
number are suppressed by inverse powers of the supersymmetry-breaking scale, which 
is thought to be at most 1 TeV. In particular, in the presence of fl-e {Vf^-Ve) niixing, the 
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diagrams of Fig. |I] are generated, which are isomorphic to the corresponding flavour- 
conserving diagrams contributing to 5a^. If the dominant supersymmetric contribution 
to (7^ — 2 comes from the chargino-sneutrino diagram Fig. |l|(b) involving left-handed 
leptons, one expects the chargino diagrams also to dominate /z — 67 H. Taking the 



BNL E821 measurement of the muon anomalous magnetic moment [12| at its face 
value fixes the overall mass scale of the sparticles circulating in the loops in Fig. |I|, and 
a supersymmetric GUT model of v^-Ve mixing can be used to calculate the amount 
of /2-e and v^-Ve mixing, i.e., the factor e in (^), enabhng the rate for — > 67 to be 
predicted. 




^ e 



5^(n) 

(a) 



(c) 



(b) 




Figure 1: Generic Feynman diagrams for fi ^ cy decay: I represents a charged slepton 
(a) or sneutrino (b), andx^"'^ and x^'^^ represent neutralinos and charginos respectively. 

Connections to other processes violating charged-lepton flavour, such as /i — > 3e and 
/i — > e conversion on nuclei, can be made in a similar way, though less directly. Among 
these processes, /i — >^ e conversion is particularly promising, since the present experimen- 
tal sensitivity may be improved by several orders of magnitude in future experiments 
such as MECO and PRISM f^. These processes receive important contributions 



from photon and Z 'penguin' diagrams, which are related to those for /i — »• 67 via a 
virtual gauge boson coupling to an e+e~ or a quark-antiquark pair, but also from box 
diagrams and their supersymmetric analogues. The dominant photonic contribution 
yields B{fiTi cTi) ^ 5.6 x 10^'^ B{fi — > 67), but this ratio may receive substantial 
corrections from subdominant contributions, which we take into account in our calcula- 
tions. Finally, we also mention the possibility in supersymmetric theories of observing 
charged-lepton- flavour violation in /xe decay [^. This involves the mixing of 

squarks as well as that of sleptons, thus providing additional information. The rate 
is small in the CMSSM with right-handed neutrinos, but might be observable using a 
future intense proton source, if tan /3 is large f\. 

In all these processes, the magnitudes of the rates depend on the masses and mixings 
of sparticles. Excessive rates for charged-lepton-flavour violation are generically pre- 



* When the right-handed sleptons also contribute significantly to the lepton-flavour-violating masses, 
the neutralino-slepton diagrams of Fig. |^(a) also contribute to — s- 67, modifying the correlation 
between ^ ^ ej and g^ — 2. 

^Larger rates are also possible in supersymmetric models with broken R parity a possibility not 
considered in this paper. 
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dieted in models with non-universal scalar masses at the GUT scale. Thus we consider 
constrained MSSM (CMSSM) models that respect this universality, such as minimal 
supergravity (mSUGRA) [|l^ , gauge- mediated supersymmetry [|l9l and no-scale mod- 



els I^Oj. In such models, even though there are no off-diagonal contributions to the 
sfermion mass matrix at Mqut-, renormalization effects due to lepton Dirac Yukawa 
couplings within the Minimal Supersymmetric Standard Model (MSSM) with massive 
neutrinos spoil this diagonal form ||2l]], making these processes observable. 



In this paper, we calculate these quantum corrections in the context of the most natural 



mechanism for obtaining sub-eV neutrino masses, namely the see-saw mechanism 
This involves Dirac neutrino masses of the same order as the charged-lepton and 
quark masses, and heavy Majorana masses M^^, leading to a light effective neutrino 
mass matrix: 

rueff = my^- {My^)-^ -ml^. (5) 

Neutrino-flavour mixing may then occur through either the Dirac and/or the Majorana 
mass matrices, which may also feed flavour violation through to the charged leptons. 

In general, the Dirac neutrino and charged-lepton Yukawa couplings, X^^^ and respec- 
tively, cannot be diagonalized simultaneously. Since both these sets of lepton Yukawa 
couplings appear in the renormalization-group equations, the lepton Yukawa matrices 
and the slepton mass matrices cannot be diagonalized simultaneously at low energies, 
either. In the basis where and are diagonal, the slepton-mass matrix acquires 
non-diagonal contributions from renormalization at scales below Mqut, of the form: 

^^l °^ Y^(2 + ~§^^^«^'^«"^3/2, (6) 

where a is related to the trilinear mass parameter: Ai = amo, where tuq is the common 
assumed value of the scalar masses at the GUT scale. 

Different oscillation scenarios for the atmospheric and solar neutrino deficits [^], e.g., 
those with small/large neutrino mixing angles and with eV or much lighter neutrinos, 
predict in general different rates for lepton-flavour violation. Typically, the larger the 
Ufj^-Ue mixing and the larger the neutrino mass scale, the larger the rates. Thus, models 
of degenerate neutrinos with bimaximal mixing lead to significantly larger effects than, 
for instance, hierarchical neutrinos with a small vacuum mixing angle. Just-so vacuum 
solutions to the solar neutrino deficit with Srn^ ~ 10~^° eV^ typically predict small 
rates if the neutrino masses are hierarchical, even if the (1-2) mixing angle is large. 



2 Sample Models of Neutrino Masses and Mixing 

In order to illustrate our estimates of the expected effects, we calculate the rates for 
rare processes violating charged-lepton number in representatives of two different types 
of models, one with small and one with large fi-e mixing. The first model (A) is 
based on Abelian flavour symmetries and symmetric fermion mass matrices [Q, and 
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leads to the following pattern of cliarged-lepton masses m^, neutrino Dirac masses mj,^, 
charged- lepton mixing Vi and Dirac mixing V^^ p5[| : 



rrii oc 
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(8) 



where e is a (small) expansion parameter related to the Abelian symmetry-breaking 
scale. In this model, even a charged-lepton matrix with large (2-3) mixing always 
predicts small /x — e mixing, as a result of fixing the charged-lepton mass hierarchies 

As a second possible model (B), we discuss the case of bimaximal mixing appearing 
In the basis where the charged-lepton mass matrix is diagonal, the neutrino 



m 



26 



mixing matrix is: 



2 
1 

V 2 



1 

2 
1 
2 



1 

V2 

1 

corresponding to a neutrino mass matrix of the form 





(9) 



eff 



m 




(10) 



where the mass parameters in this texture are related to the mass eigenvalues by 

m = 712^/2, 
es = (m2 + mi)/4m, 
6 = V2 {nil — 7712) / 4:771 . 



(11) 
(12) 

(13) 



In our analysis we assume hierarchical neutrinos, but the case of degenerate neutri- 
nos can be treated similarly. This would lead to different predictions, since the neu- 
trino mass scales change, typically with larger rates for charged-lepton-flavour violation. 
Hence our results for model (B) are quite conservative. 



3 Supersymmetric Calculations 

We have calculated the rates for processes violating charged lepton number in both 
frameworks, including complete sets of one-loop sparticle diagrams. We parametrize 
the universal soft supersymmetry-breaking masses by the GUT-scale parameters mo 
and 7171/2, for sfermions and gauginos respectively, and use the renormalization-group 
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equations of the CMSSM to calculate the low-energy sparticle masses pT|. Other rel 



evant free parameters of the MSSM are the trilinear coupling A, the sign of the Higgs 
mixing parameter /x, and the ratio of Higgs vev's, tan/5. Here we consider only /i > 0, 
since this is the sign favoured by (yf^ — 2. The physical charged-slepton masses are found 
by numerical diagonalization of the following matrix: 



where all the entries are 3x3 matrices in flavour space. Using the superfield basis 
where the Yukawa coupling matrix is diagonal, we can write: 

^II = {mjy + Sml^+mj- Ml{^- sin'^ew)cos2p (15) 

= (mf-^)2 + m^-M|sin2^VKCOs2/5 (16) 

^RL = (^f + Me-/xtan/5)m^ (17) 

mln = rn^^ (18) 

where (m|)^, (^e^)^ denote the diagonal contributions to the corresponding 

matrices, obtained by numerical integration of the renormalization-group equations, 
and ^rr?y^ and bAi denote the off-diagonal terms that appear because and \i may 
not be diagonalized simultaneously - see (^. Analogously, for sneutrinos we have 

= {mlf + Sml^ + ^Mlcos2P (19) 

The mixing parameter e in (^ is given in terms of the parameters of this matrix. In 
the simplified case where the lepton-number-violating mass terms are in the iL sector, 
so that the chargino-sneutrino diagram dominates both — 2 and fi 67, assuming 
that the sparticle masses have (approximately) a common value m, and making a naive 
mass-insertion approximation, one would find 



(m^ 

/m 2 



(20) 



but this is only indicative, and we use complete formulae in our results below. 



We start by fixing the elements of the Yukawa coupling matrices at the GUT scale to be 
consistent with the experimental values of the fermion masses and the absolute values 
of the CKM matrix elements 0. This is done by choosing appropriate coefficients of 
order one in the entries of the lepton matrices. In the notation of 0], we choose for 
model (A): 

Cu = 0.77, C23 = 0.79, (21) 

and for model (B): 

Cu = 2.75, C23 = 1.13, (22) 

with the unspecified coefficients taken as unity. These coefficients do not change signif- 
icantly for the two values of tan j3 considered in the present work. 
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We then use the full Yukawa coupling matrices in the renormalization-group equations 
including the effects of X^^ on the CMMSM parameters at the see-saw mass scale 



We have checked that our results using the full matricial forms for the Yukawa couplings 
do not differ significantly from the common approach of considering diagonal Yukawa 
matrices and neglecting the lighter generations. We check our results by constructing 
the slepton mass matrices in the superfield base where the Xi are diagonal and inserting 
the non-diagonal elements induced by the presence of A^^ on the renormalization-group 
equations between the GUT and see-saw mass scales. Finally, we use the full matricial 
forms for all the parameters which appear in the vertices in the diagrams of Fig. |l|, and 
the results of |28| to calculate the rates for /i — 67 and /i — >■ e conversion, and those 
of P to calculate the rate for Kl — >• /le. 



4 Constraints on the CMSSM 

We display in the remaining figures the {mQ,mi/2) planes in the CMSSM for tan/? = 
10, 30, assuming yU > as suggested by the sign of 6a^, and Aq = 0. The experimental 
constraints on the CMSSM are taken from |^9|, |13|, where further details of their imple- 
mentation can be found. We note in particular that the following choices are used here 
for the pole mass of the top quark: rrit = 175 GeV, and for the running mass of the bot- 
tom quark: mh{mi,)^^ = 4.25 GeV. We combine the constraints given in |TB[ with 
the contours suggested by the neutrino mass textures introduced above for B{fi 67) 
in Fig. |, for Biix'Ti e'Ti) in Fig. |, and for B{Kl fx^e^) in Fig. § 

The dark (brick-red) shaded regions in the {mo,mi/2) planes in these figures are ex- 
cluded 1^ because the lightest supersymmetric particle is the lighter f, which is dis- 
allowed by the astrophysical requirement that cold dark matter be electrically neutral. 
The light (turquoise) shaded regions are those where the LSP is the lightest neutralino 
and its cosmological relic density f^^^^ ^^^^ favoured range 0.1 < Q^h"^ < 0.3 



Lower values of f^^^^ would be possible if there are other sources of cold dark matter, 
whereas larger values of Q^h"^, which occur generically at larger values of mi/2 and mo, 
are excluded by cosmology. 

We display as (red) dash-dotted lines mass contours for the lightest CMSSM Higgs 



boson: nih = 113,117 GeV, as calculated in ||2^. This range corresponds roughly 
to values of mi/2 and mo that are compatible, within theoretical errors, with the LEP 
Higgs 'signal' at uih = 115^1'^ GeV ||30[, for our default choices of Aq, rut and mb{mb)^M- 
There is good overall consistency between mh and the other constraints for 10 < tan j3 < 
55, but we do not insist on the range 113 GeV < rrih < 117 GeV, in view of the 
theoretical uncertainties and because the LEP Higgs 'signal' might turn out to be a 
false alarm, in which case rrih could be larger. 

The medium dark (green) shaded regions in panels (b) and (d) are excluded by our 
implementation of the 6 — > 57 constraint. As described in we use the latest NLO 
QCD calculations for large tan/5 and allow values of mi/2 and mo that, after including 
the expected theoretical errors due to the scale choice and model dependences, may fall 
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within the 95% confidence level range 2.33 x 10"^ < B{b ^ s < 7) < 4.15 x 10"''. In 
the panels (a) and (c) of the figures, for tan/3 = 10, there is no relevant constraint from 
6 — > 57, and we display as a dashed line the LEP lower limit m^± > 104 GeV. There 
is a similar constraint in the panels (b) and (d), for tan/? = 30, which is omitted for 
clarity. 



To complement our summary of [2£], we note that the regions allowed by the E821 



measurement of at the 2-a level [|T2l are shown in the figures as light (pink) shaded 



regions with solid black line boundaries |13]. Also shown as black dashed lines are 
the regions favoured by at the 1-a level. We emphasize the impressive consistency 
between the constraints from a„, rrih, b ^ S'j and cosmology for tan/? > 10. As 



discussed in [0, combining all the other constraints with the 1-a range for a^, one 
finds quite small allowed regions of the (mi/2, mo) plane centred on: ~ (250, 100) GeV 
for tan/? = 10 (see panels (a) and (c)) and ~ (350, 170) GeV for tan/3 = 30 (see panels 
(b) and (d)). 



5 Results for Processes Violating Charged-Lepton 
Number 



Fig. 1^ displays the predictions for B{fi — >• 67) in texture (A) (panels (a) and (b)) 
and texture (B) (panels (c) and (d)). We see that texture (A) generically predicts 
that B{fi —>■ 67) should occur within one or two (two or three) orders of magnitude 
of the present experimental upper limit if lies within the l(2)-cr range suggested by 
E821 ll^. Within this model, the experimental upper limit on B{fi cj) < 1.2 x 10~^^ 
excludes a domain of the (mi/2, mo) plane, close to the origin, that may be compared to 
that excluded by slepton searches at LEP. When tan (3 ~ 10, it also has a similar effect 
to that of the upper limit on the supersymmetric contribution to — 2. For larger 
tan P ~ 30, the constraint due to the present upper limit on B{fi 67) is intermediate 
between the g^ — 2 and b ^ s'-y constraints. In model (B), we find values of B{fi 67) 
that are characteristically about an order of magnitude smaller than in model (A) in 



the parameter region allowed by cosmology ||29| and g^ — 2 [|T^]. 



As already mentioned, the rate for ^~Ti —>■ e^Ti conversion is linked to that for 
B{fi — ^ 67), with a proportionality coefficient ~ 5.6 x 10~^ if fi~Ti — > e~Ti is dominated 
by photon exchange. However, this is not necessarily the case, since Z° exchange and 
box diagrams may also contribute, rendering the ratio B{^~Ti — > e^Ti)/B{ji — > 67) 
non-universal. Accordingly, we plot in Fig. ^ the predictions of texture (A) (panels (a) 
and (b)) and texture (B) (panels (c) and (d) for B{fi~Ti — > e~Ti). We see that there 
are large domains of the (mi/2,mo) plane where these textures suggest that B{fi~Ti 
e~Ti) > 10~^^, which is the sensitivity of the proposed MECO experiment |jl^, |^. In 
model (A), these include all the regions allowed by g^ — 2 at the 2-a level, and, in model 
(B), most of the allowed region. We infer that the physics interest of this proposal is 
greatly enhanced by the recent result on g^ — 2 from E821 [T^. 
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m,,2 (GeV) mi/2 (GeV) 



Figure 2: The contours B{ii — *• 67) = 10"^"^, 10~^^, lO^"*^^ and lO""*^^ in (a, b) texture 
(A) and (c, d) texture (B) are shown as dash-dotted black lines in the (mi/2,mo) planes 
for /i > and tan/5 = (a, c) 10 and (b, d) 30. Other constraints in these planes are 
taken from assuming Aq = 0,mt = 175 GeV and mb{mb)^M = 4.25 GeV. The 
regions allowed by the E821 measurement of at the 2-a level / |7^/ are shaded (pink) 
and bounded by solid black lines, with dashed lines indicating the 1-a ranges jl^/ . The 
dark (brick-red) shaded regions are excluded because the LSP is the charged fi, and the 
light (turquoise) shaded regions are those with 0.1 < Vl^h'^ < 0.3 that are preferred by 
cosmology. We show the contours rrih = 113, 117 GeV, and in panels (a, c) the contour 
m^± = 104 GeV. The medium (dark green) shaded regions are excluded by b ^ s'j. 
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mi/2 (GeV) mi/2 (GeV) 



Figure 3: The contours Bifi'Ti e^Ti) = IQ-^^, 10"^^, 10-^^ lO^^^ and lO^^^ m 
(^a, b) texture (A) and (c, d) texture (B) are shown as dash-dotted black lines in the 
{mi/2,mo) planes for fi > and tan/? = (a, c) 10 and (b, d) 30. Other constraints in 
these planes are taken from 



/, as described in the caption of Fig. 
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Finally, we plot in Fig. ^the corresponding model predictions for B{K^ /i^e^). This 
process is very interesting |p, because it combines flavour violation in the quark and 
lepton sectors. For this same reason, one expects rather small values of B{Kl ^^e^), 
far below the present experimental upper limit. However, again for this same reason, 
there is clearly even more uncertainty in the predictions for B{K^ /^^e''') than there 
was already for B{fi — >• ej) and B{fj,~Ti e~Ti). The sensitivity to the lepton mass 
texture is seen clearly by comparing panels (a, b) and (c, d) of Fig. ^. Nevertheless, 
we note that B{K^ —>■ yU^e^) > 10^^^ in a significant fraction of the parameter region 
favoured by (7^ — 2 at the one-cr level. Thus, we think that this process is of potential 
interest at an intense proton source. 



6 Conclusions 



We have argued in this paper that the BNL E821 measurement of — 2 |]12|, taken 



at face value, may be used to normalize predictions for the charged-lepton-number- 
violating processes /i cy, /i — e conversion and /^^e^, within a supersymmetric 

GUT framework. We have illustrated our argument with a couple of specific textures 
for fermion masses that are consistent with the data on neutrino oscillations. In these 
examples, we find that fi ^ cy decay may appear at a rate within one or two (two or 
three) orders of magnitude of the present experimental upper limit if — 2 lies within 
its present one- (two-) a range. These models also make us optimistic that fi ^ e 



conversion on heavy nuclei may be accessible to the next round of experiments [^, |^ . 
The prospects for observing — > fi'^e^ decay are not so rosy, but this decay might 
also be accessible to some future round of experiments with an intense proton source. 

The discovery of neutrino oscillations has been a major breakthrough in flavour physics. 
If confirmed, the deviation oi — 2 from the SM prediction would be a breakthrough 
towards new physics at the TeV scale. Their combination suggests not only that the 
conservation of charged lepton number is not sacred, but also that its violation may 
soon be observable. If so, this would be an invaluable new window on the physics 
of lepton fiavour, as well as on physics at the TeV scale. It would provide a bridge 
between neutrino oscillations and accelerator physics, as well as yield novel information 
on lepton mixing. 

We encourage strongly the most sensitive possible experiments to probe the violation of 
charged lepton number. 
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Figure 4: The contours B{Kl /i±e^) = 10-l^ 10-^°, lO^^i and lO^^s ^a, 5^ 
texture (A) and (c, d) texture (B) are shown as dash-dotted black lines in the {mi/2,mo) 
planes for /i > and tan/? = (a, c) 10 and (h, d) 30. Other constraints in these planes 
are taken from J^, [7^/ , as described in the text. 
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discussions related to this analysis. 
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